Modeling long-range memory with stationary Markovian processes 
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In this paper we give explicit examples of power-law correlated stationary Markovian processes y(t) 
where the stationary pdf shows tails which are gaussian or exponential. These processes are obtained 
by simply performing a coordinate transformation of a specific power-law correlated additive process 
x(t), already known in the literature, whose pdf shows power-law tails x~ a . We give analytical and 
numerical evidence that although the new processes (i) are Markovian and (ii) have gaussian or 
exponential tails their autocorrelation function still shows a power-law decay (y(t)y(t + r)) oc r - * 3 
where (3 grows with a with a law which is compatible with /3 = a/2 — rj. When a < 2(1 + rf) the 
process y(t), although Markovian, is long-range correlated. Our results help in clarifying that even 
in the context of Markovian processes long-range dependencies are not necessarily associated to the 
occurrence of extreme events. Moreover, our results can be relevant in the modeling of complex 
systems with long memory. In fact, we provide simple processes associated to Langevin equations 
thus showing that long-memory effects can be modeled in the context of continuous time stationary 
Markovian processes. 

PACS numbers: 02.50.Ey, 05.10.Gg, 05.40.-a, 02.50.Ga 



I. INTRODUCTION 

Stochastic processes have been used to model a great 
variety of systems in disciplines as disparate as physics [l|, 

f, m IE E| , genomics 0Jp , finance [1, , climatology 
and social sciences [12 ]. 

One possible classification of stochastic processes takes 
into account the properties of their conditional probabil- 
ity densities. In this respect, Markov processes play a 
central role in the modeling of natural phenomena. In 
the framework of discrete time stochastic processes, a 
process x(t) is said to be a Markov process if the condi- 
tional probability density P(x n +i, t n +i\x n , t n ; . . . ]X\,t\) 
depends only on the last value x n at t n and not on the 
previous values x n -i at t n -\, x„_2 at t n -i, etc. More 
generally, the transition probability of any Markov pro- 
cess fulfills the Chapman-Kolmogorov equation [2(]. It 
is worth noting that a Markov process is fully deter- 
mined by the knowledge of the probability density func- 
tion (pdf) W(x, t) of the process and the transition prob- 
ability P(x n+ i,t n+ i\x n ,t n ). When the Markovian pro- 
cess is continuous both in space and time, the time evo- 
lution of the pdf is described by a Fokker-Planck (FP) 
equation. Such level of simplicity is rather unique among 
stochastic processes. In fact, a non-Markovian process is 
characterized by an infinite hierarchy of transition prob- 
abilities. In this case, the time evolution of the pdf is 
described by a Master Equation rather than a simpler 
FP equation. 

Another classification of stochastic processes considers 
the nature of correlation of the random variable. Under 
this classification, random variables arc divided in short- 
range and long-range correlated variables. Short-range 
correlated variables are characterized by a finite mean of 
time-scales of the process whereas a similar mean time- 
scale does not exist for long-range correlated variables 



[131 ]. An equivalent definition can be given by consider- 
ing the finiteness or infiniteness of the integral of the au- 
tocorrelation function of the random process [3, [H, EH • 
In the presence of long-range correlation, the time inte- 
gral s(t) of the process x(t) is a super-diffusive stochastic 
process showing (|As(t)| 2 ) ~ D 1 i 7 where 7 > 1 and D 1 
is a constant. Superdiffusive stochastic processes have 
been observed in several physical systems. A classical 
example is Richardson's observation that the relative sep- 
aration £ of two particles moving in a turbulent fluid at 
time t follows the relation < £ 2 (t) >cx t 3 [13 • Other 
examples include anomalous kinetic in chaotic dynamics 
due to flights and trappi ng Jl8l| . dynamics of aggregate 
of amphiphilic molecules [19[ , dynamics of a tracer in a 
two-dimensional rotating flow [20( , non-coding regions of 
complete genomes [2lj and volatility in financial markets 

Several stationary Markovian processes are short-range 
correlated. In fact, the paradigmatic Markovian process 
is the Ornstein-Uhlembeck (OU) one [23j], whose auto- 
correlation function is the exponential function p(r) = 
e _T / T where T is the time-scale of the process. Although 
in the OU process there is one single time-scale, a general 
Markovian stationary process can be multi-scale, i.e. it 
may admit either a discrete or a continuum set of time- 
scales. In the last case, when the largest time-scale is 
removed to infinity the process can even be long-range 
correlated. The paradigmatic Markovian process with 
power-law autocorrelation function is given by the fam- 
ily of processes considered in Ref. [2J|. These are sta- 
tionary Markovian power-law correlated processes that 
were introduced in the context of diffusion in optical lat- 
tices and semiclassically describe the motion of atoms in 
a one-dimensional optical lattice formed by two counter- 
propagating laser beams perpendicularly polarized. For 
a certain choice of the relevant paramenters the processes 



2 



become long-range correlated. 

The existence of a power-law decaying autocorrelation 
function in the processes of Ref. [24[ is intimately related 
to the existence of power-law tails in the stationary pdf. 
This is easily understood by considering that the pro- 
cesses of Ref. [24j describe particles moving in a confin- 
ing Smoluchowski potential which asymptotically grows 
like log(x). If one compare such slow growth with the 
one associated to the OU process, whose Smoluchowski 
potential grows like x 2 , it is easy to recognize that in 
the case of Ref. {24| (i) a particle can reach positions far 
away from the center of the potential because it is subject 
to a relatively weaker force and (ii) if a particle reaches 
a position X, then it is not suddenly recalled towards 
the center of the potential and therefore it can explore 
for relatively long times the regions around X. Loosely 
speaking, the time-series of the processes of Ref. [24| can 
show persistencies and clustering of extreme events. Such 
processes perfectly fit the features of the model proposed 
in Ref. [25|], where long-range dependencies are shown 
to explain the clustering of extreme events. However, 
one could have in principle slowly decaying autocorrela- 
tion functions without necessarily observing the occur- 
rence of extreme events. One such example is given by 
the Fractional Brownian motion (FBm) [13| which is a 
stochastic process where the autocorrelation function de- 
cays like a power-law and the stationary pdf is gaussian. 
In this paper, in the context of Markovian processes, we 
give explicit examples of power-law correlated stationary 
processes where the stationary pdf shows tails which are 
gaussian or exponential. We will introduce such processes 
starting from appropriate coordinate transformations of 
an additive processes introduced in Ref. [2J|. 

The paper is organized as follow. In section |TT] we re- 
view the eigenfunction methodology used to analyze the 
correlation properties of a given stochastic process and 
introduce a specific power-law correlated process with 
power-law tails. In section Unl and UVl we will present ex- 
amples of power-law correlated stochastic processes with 
gaussian and exponentail tails in the stationary pdf re- 
spectively. In section [V] we will draw our conclusions. 



II. POWER-LAW TAILS IN THE PDF 

In this section we will briefly review the family of 
stochastic processes introduced in Ref. [2J] and whose 
ergodicity properties have been investigated in Ref. [2(| . 
A similar class of such processes have been considered in 
Ref. [13. 

Let us consider a continuous Markovian stochastic 
process x(t) whose pdf W(x,t) is described by the 
FP equation with constant diffusion coefficient dtW = 
-d x {D^\x)W) +Dd 2 x W. For the sake of simplicity, in 
this study we set D = 1. In general, the eigenvalue spec- 
trum of the FP equation describing a stationary process 
consists of a discrete part Ao = 0, Ai, X p and a con- 
tinuous part ]A c ,+oo[ (A c > X p ) associated with eigen- 



functions (p\. The stationary pdf is W(x) = ipo. The 
FP equation with constant diffusion coefficient can be 
transformed into a Schrodinger equation [3[ with a quan- 
tum potential V s (x) = (D& (x)) 2 /4 + d x D^ (x)/2. The 
eigenvalue spectrum of the Schrodinger equation is equal 
to the eigenvalue spectrum of the FP equation. The 
relation between the eigenfunctions of the FP equation 
and the eigenfunctions ip\ of the Schrodinger equation is 
fx = V-'aV-'o- For a stationary process the 2-point proba- 
bility density function W 2 (x, t\x' ,t + r) can be expressed 
in terms of the eigenfunctions of the Schrodinger equa- 
tion. Specifically, one can write 



W 2 (x,t;x',t + T) = tp {x) ipo(x') x 

X P P + OC 

,A=Ai ^ A <= 



(1) 



-At 



Eq. {T]) extends the analogous expression valid for a FP 
equation with only discrete spectrum [3[ to the case in 
which there also exists a continuous part of the spectrum. 
By direct inspection, it can be shown that W 2 fulfills 
the Chapman-Kolmogorov equation. In order to evalu- 
ate the autocorrelation function p(r) — ((x(t + r)x(t)) — 
(x(t)) 2 )/({x 2 (t)) - (x(t)) 2 ) of the stochastic variable x(t), 
we make use of the expression 



(x(t + r)x{t)) 



A=Ai 



C 2 x e- XT d\, 



(2) 



where C\ = J dx x(p\(x). Eq. ^ follows from Eq. (fT]) 
and from the definition of the autocovariance function 



(X(t + T)x(t)) 



dx' dx x' xW 2 (x,t 1 x' ,t + t). (3) 



Eq. @ holds true under the assumption that the inte- 
grations in J dx' J dx and J dX can be interchanged. 

The asymptotic temporal dependence of the auto- 
correlation function can have a different behavior con- 
ditioned by the properties of the eigenvalue spectrum 
[H, [2^, [lo, HH . Specifically, following [32[ one can dis- 
tinguish three different cases, depending on the existence 
of a continuum spectrum of eigenvalues and whether or 
not such spectrum is attached to the ground state. 

In fact, the class of processes introduced in Ref. [24| 
belongs to the one admitting a continuum part of the 
spectrum attached to the ground state. In this paper we 
will consider the specific stationary Markovian processes 
associated with a quantum potential Vs given by 



V s (x) = 



-V if \x\ < L, 
Vi/x 2 if |a;| > L, 



(4) 



where L, Vo and V\ are positive constants. The reason 
for considering such specific potential, among all those 
fullfilling the requirements of Ref. [24j ]. is that it is ex- 
actly solvable and therefore it will allow us to perform 
most calculations analytically. 



The parameters L, Vq and V\ can be chosen in such a 
way that the spectrum contains one single discrete eigen- 
value Ao = and a continuous part for A > 0. As a result, 
the parameters L, Vq and Vy are not independent. In 
fact, the continuity of d x tpo in x = L provides a relation 
between them. The Langevin equation of the process is 



x = h(x) + dz 

( -2yJ% tan(/^a;) if |x| < L, 

Hx) = { 

[ (1 - VTTTV[)/x if |x| > L. 

V x = J Ltan( v / Vo"i)fl + itan(VVo"i) 



(5) 



The associated FP equation describes the dynamics of 
an overdamped particle moving in a Smoluchowski po- 
tential U[x) — — J dx h(x) that increases logarithmically 
in x. For |x| ^ L, the eigenfunction of the ground state 
is tpo = Bq cos(y / V u ":r) whereas for |x| > L it decays ac- 

r (l-Vl+4Vi)/2 



cording to ip — A 



The constants A$ and 



B are set by imposing that ^0 is normalized and con- 
tinuous in x — L. It is worth noting that for |x| > L the 
stationary pdf W(x) of the stochastic process is a power- 
law function decaying as \x\~ a with a = VI + 4 V\ — 1. 
The normalizability of the eigenfunction of the ground 
state is ensured if a > 1. In the present study we con- 
sider stochastic processes with finite variance which im- 
plies a > 3. Due to parity arguments, only the odd 
eigenfunctions %j:^ dd ^ of the continuous spectrum give 
a non- vanishing contribution to C\. For |x| > L the 

eigenfunction ip^ dd ^ is a linear combination of Bessel 
functions t()^ dd ' — A\yfxJ v (yf\x 



B X y/xY v {V\x) 



where v 



(a 



l)/2. For \x 
D\ sm(y/Vo + As). The coefficients A\ 

1 (odd) 



^ L we find ip^ dd ^ 



B\ and D\ are 

fixed by imposing that ip\^ aU/ ' and its first derivative are 

continuous in x — L and that il)^ dd ^ are orthonormalized 
with a (5-function of the energy. Similar conditions apply 
to the even solutions. 

By using these eigenfunctions we obtain an exact ex- 
pression for C\. The further integration required in Eq. 
([2]) to obtain (x(t + r) x(t)) cannot be performed analyt- 
ically. By using Watson's lemma [33| and by considering 
that the first term of the Taylor expansion of C\ is pro- 
portional to A^" -5 )/ 2 for small values of A, for large values 
of r one gets (x(t + r) x(t)) cx t~P where (3 = (a — 3)/2. 
That indicates that this stochastic process is stationary, 
Markovian and asymptotically power-law autocorrelated. 
When 3 < a < 4 the process is long-range correlated. 



III. GAUSSIAN TAILS IN THE PDF 

In this section we explicitely present a stationary 
Markovian process with power-law decaying autocorre- 
lation function and a stationary pdf with gaussian tails. 




FIG. 1: The Figure shows the drift coefficient H(y) (top 
panel) and diffusion coefficient G(y) (bottom panel) of the 
process defined by the coordinate transformation of Eq. (|6]) 
for the case when L = 1.0, Vo = 0.987 (i.e. a = 3.05) and 
s = 1.0 (i.e. r = 1.2096). 



In fact, let us consider the coordinate transformation: 



12s Erf 



(1 - r(x))Erf(£r/V2s)+ 



1M 



a(x)(l + v) 



\x\ > L 



Id < L 



L r y/Va sec(L y/Vo) 2 + r tan(L y/Vo) 



r(x) 



2A\(Lx a - xL a ) 
(a — 1) L a 1° 



(6) 



where r is a real positive constant. By using the Ito 
lemma, one can show that, starting from the process of 
Eq. ([5]), in the coordinate space y = f g (x) one gets a 
multiplicative stochastic process whose stationary pdf is: 



Nt e 



w g {y) 



\y\> Lr 



N n cos(\/Voy/r) \y\< Lr 



(7) 



where iVj and Nn are normalization constants that can 
be analytically computed by imposing that W g (y) is con- 
tinuous in y — ±Lr and it is normalized to unity. The 
real constant r is fixed by imposing that the diffusion co- 
efficient G(y) of the multiplicative stochastic process in 
the y coordinate space is continuous in y = ±Lr. In Fig. 
Q]we show the drift coefficient H(y) (top panel) and dif- 
fusion coefficient G(y) (bottom panel) for the case when 
L = 1.0, Vo = 0.987 (i.e. a = 3.05) and s = 1.0 (i.e. 
r = 1.2096). The diffusion coefficient G(y) is continuous 
in y = ±Lr although its first derivative is discontinu- 
ous. The drift coefficient H(y) suffers a discontinuity in 
y = ±Lr. 
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FIG. 2: In the top panel we report the results of the numerical 
integration of Eq. J8J for the case when L = 1.0, s — 1.0 and 
the Vb values are choosen in such a way that the parameter a 
assumes the values shown in the legend. In the bottom panel 
we report the values of the exponents f3 g obtained by per- 
forming a nonlinear fit of the autocorrelation function shown 
in the top panel. 



formed for the case when L — 1.0, Vb = 0.987 (i.e. 
a = 3.05) and s = 1.0 (i.e. r = 1.2096). The solid 
(red) line shows the theoretical prediction obtained from 
Eq. ([5]), while the open circles show the result of the 
numerical simulations. By performing a nonlinear fit 
(dashed blue line), the autocorrelation function shows 
an asymptotic decay compatible with a power-law T~P g , 
with (3 g = 0.86. In the inset of the top panel we show 
the numerical simulation (circles) relative to the mean 
square displacement (\As(t)\ 2 ) where s(t) is the stochas- 
tic process obtained by integrating over time the process 
defined by the coordinate transformation of Eq. ([5]). A 
nonlinear fit (solid blue line) shows that (|As(t)| 2 ) oc t s 
with 5 = 1.21, thus confirming that we are observing 
a superdiffusive long-range correlated stochastic process. 
The bottom panel of Fig. [3] shows the stationary pdf of 
the process. Again the solid (red) line shows the theo- 
retical prediction of Eq. ([7]), while the open circles show 
the result of the numerical simulations. In the inset we 
show the same pdf in a shorter range of values in order 
to emphasize that inside the region |y| < L r the pdf has 
a behaviour different from gaussian. 

Those shown in Fig. [3] are time-average numerical sim- 
ulations performed according to the relation: 



The autocorrelation function of the process defined by 
the coordinate transformation of Eq. © is given by 
p g (r) = ({y{t+T)y{t))-{y{t))*)/{{y\t))-{y{t))Z) where: 

(y(t)y(t + r)) = / d\C 2 x e- XT 
Jo 

(8) 

/+oo 
dx f g (x)tpo(x) ip\(x) 
-oo 

where ipo(%) and ip\( x ) are t ne eigenfunctions of the pro- 
cess of Eq. ([5]). Eq. can be used to numerically 
obtain the autocorrelation of the process defined by the 
coordinate transformation of Eq. ^ . 

In the top panel of Fig. [5] we report the results of 
the numerical integration of Eq. ([5]) for the case when 
L = 1.0, s = 1.0 and the Vb values are choosen in such 
a way that the parameter a assumes the values shown 
in the legend. The asymptotic behaviour of these auto- 
correlation functions seems compatible with a power-law 
r - ^ 9 . In the bottom panel of Fig. [2] we report the values 
of the exponents p g obtained by performing a nonlinear 
fit of the autocorrelation function shown in the top panel. 
Such values show a dependence from the a parameter 
which seems compatible with a linear law [3 g = a/2 — r/ g , 
with r\ g ss 0.61. For values a < 2-q g + 1 we get (i g < 1, 
i.e. the stochastic process thus generated is long-range 
correlated. It is worth noticing that the autocorrelation 
function (y(t)y(t + r)) does not show any dependance 
from the s parameter. 

In the top panel of Fig. [3J we show the results of nu- 
merical simulations of the autocorrelation function per- 



1 r 

P {t) = - dt x*(t)x*(t + r) (9) 
1 Jo 

where T is the length of the simulated time-series and 
x*(t) is one realization of the process. Indeed, in order 
to improve the statistical reliability of our numerical sim- 
ulations, in the region r ^ 1 we have also averaged over 
a number M of different realizations of the process: 

1 M 1 r T 

P t ^=mY,t dtx 3 {t)x {t + T) (10) 

The data shown in the figure are the mean and the stan- 
dard deviations of the M autocorrelation values com- 
puted in each iteration for each time lag. The values 
of M are M = 10 in the region r e [1, 10] and M = 20 
in the region r ^ 10. The size of each time-series was 
T = 1.01 10 11 with a time-step of At = 0.005. The start- 
ing points of the simulated time-series were all the same 
with Xj(0) = 0.1 where j = 1, • • • , M. In order to sim- 
ulate the process in the y coordinate space, we start by 
simulating the process of Eq. ([5]) and compute y — f(x) 
for each simulated x value. However, we have explicitely 
checked that such procedure is equivalent to a direct sim- 
ulation of the Langevin equation obtained starting from 
H(y) and G{y). 

The existence of power-law correlated processes with 
gaussian tails does not contraddict the Doob Theorem 
j34J | . In fact, such theorem deals with the case when the 
process admits stationary pdfs and 2-point conditional 
transition probabilities which are both non singular and 
gaussian on the whole real axis. 
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by: 



FIG. 3: The figure shows time-average numerical simulations 
performed according to Eq. (|10[) performed for the case when 
L = 1.0, Vo = 0.987 (i.e. a = 3.05) and s = 1.0 (i.e. 
r — 1.2096). The simulation parameters are: M = 10 in 
the region r £ [1, 10] and M = 20 in the region r ^ 10, 
T = 1.01 10". The time-step was At = 0.005. In the top 
panel we show the results for the autocorrelation function. 
The solid (red) line shows the theoretical prediction obtained 
from Eq. (|8]), while the open circles show the result of the 
numerical simulations. By performing a nonlinear fit (dashed 
blue line), the autocorrelation function shows an asymptotic 
decay compatible with a power-law r~^ a , with /3 g = 0.86. In 
the inset of the top panel we show the numerical simulation 
(circles) relative to the mean square displacement (|As(i)| 2 ). 
A nonlinear fit (solid blue line) shows that (|As(i)| 2 ) oc t s 
with 8 = 1.21. The bottom panel shows the stationary pdf of 
the process. Again the solid (red) line shows the theoretical 
prediction of Eq. (J7), while the open circles show the result 
of the numerical simulations. 



IV. EXPONENTIAL TAILS IN THE PDF 



In this section we explicitely present a stationary 
Markovian process with power-law decaying autocorre- 
lation function and a stationary pdf with exponentially 
decaying tails. In fact, let us consider the coordinate 
transformation: 



ilog(7(a:-L) + e-T £p ) \x\> L 

fe(x) = <; (11) 

rx \x\ < L 




(12) 



-7 e- 2 ^ < 1+ e t e ;i +A \y\>Lr 



2rVVotan(v%y/r) \y\<Lr 



By using the Ito lemma, one can show that, starting from 
the process of Eq. |5|), in the coordinate space y — f e (x) 
one gets the multiplicative stochastic process described 



A = 7 L - e~~< Lr 

where r is a real positive constant which is fixed by im- 
posing that the diffusion coefficient G(y) is continuous in 
y = ±Lr. It is straightforward to prove that such process 
admits the stationary pdf: 

( Ni e~™ \y\>Lr 
W e (y) = { _ (13) 

[ Nii cos(V%y/r) \y\<Lr 

whose tails are exponential. N] and Nu are normal- 
izations constants that can be analytically computed by 
imposing that W e (y) is continuous in y = ±Lr and it is 
normalized to unity. 

The autocorrelation function of the process of Eq. (fT5|) 
can be obtained starting from Eq. (JSj) with f g (x) now 
replaced by f e (x) of Eq. (fTTj) . 

In the top panel of Fig. [4] we report the results of 
the numerical integration of Eq. ||5J) for the case when 
L = 1.0, 7 = 1.0 and the Vo values are choosen in such 
a way that the parameter a assumes the values shown 
in the legend. The asymptotic behaviour of these auto- 
correlation functions seems compatible with a power-law 
t~P" . In the bottom panel of Fig. 2] we report the values 
of the exponents f3 e obtained by performing a nonlinear 
fit of the autocorrelation function shown in the top panel. 
Such values show a dependence from the a parameter 
which seems compatible with a linear law f3 e — a/2 — rj e , 
with rj e w 0.68. Differently form the gaussian case, now 
the autocorrelation function (y(t)y(t + t)) seems to show 
some dependance from the 7 parameter. As an exam- 
ple, we have computed the autocorrelation functions for 
the same values as above and with 7 = 1.0 replaced by 
7 = 10.0. Again we find that r\ depends upon a according 
to a linear law (3 e — a/2 — rj e where now r] e w 0.62. 

In the top panel of Fig. O we show the results for the 
case when L = 1.0, V = 1.020 (i.e. a = 3.21) and 
7 = 1.0 (i.e. r = 0.567). The solid (red) line shows 
the theoretical prediction obtained from Eq. ([5]), while 
the open circles show the result of the numerical simula- 
tions. By performing a nonlinear fit (dashed blue line), 
the autocorrelation function shows an asymptotic decay 
compatible with a power-law t - ^, with /3 e = 0.79. In 
the inset of the top panel we show the numerical simu- 
lation (circles) relative to the mean square displacement 
(|As(t)| 2 ) where s(t) is the stochastic process obtained 
by integrating over time the process defined by the coor- 
dinate transformation of Eq. (fTT]) . A nonlinear fit (solid 
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FIG. 4: In the top panel we report the results of the numerical 
integration of Eq. ((8} for the coordinate transformation of Eq. 

when L — 1.0, 7 = 1.0 and the Vo values are choosen in 
such a way that the parameter a assumes the values shown in 
the legend. In the bottom panel we report the values of the 
exponents /3 S obtained by performing a nonlinear fit of the 
autocorrelation function shown in the top panel. 



blue line) shows that (|As(t)| 2 ) oc t s with S = 1.26, thus 
confirming that we are observing a superdiffusive long- 
range correlated stochastic process. The bottom panel 
of Fig. [5] shows the stationary pdf of the process. Again 
the solid (red) line shows the theoretical prediction of Eq. 
(|13[) . while the open circles show the result of the numer- 
ical simulations. In the inset we show the same pdf in a 
shorter range of values in order to emphasize that inside 
the region \y\ < L r the pdf has a behaviour different 
from exponential. 

Those shown in Fig. [5] are time-average numerical 
simulations performed according to Eq. ^j. Differently 
from the previous case, when simulationg the process we 
directly consider the Langeving Equation of Eq. (TT2")) . 
Again, in order to improve the statistical reliability of 
our numerical simulations, in the region r > 1 we have 
also averaged over a number M of different realizations 
of the process, according to Eq. (fTOj) . The data shown 
in the figure are the mean and the standard deviations 
of the M autocorrelation values computed in each itera- 
tion for each time lag. The values of M are M = 10 in 
the region r 6 [1, 10] and M = 20 in the region t ^ 10. 
The size of each time-series was T = 1.025 10 11 with a 
time-step of At = 0.01. The starting points of the sim- 
ulated time-series were all the same with Xj(Q) = 0.1 
where j = 1, • • • , M. 

V. CONCLUSIONS 

In summary, we have shown new stationary Markovian 
processes which are power-law correlated and have a sta- 



FIG. 5: The figure shows time-average numerical simulations 
performed according to Eq. (| lOf) performed for the case when 
L = 1.0, Vo = 1.020 (i.e. a = 3.21) and 7 = 1.0 (i.e. 
r = 0.567). The simulation parameters are: M = 10 in 
the region r € [1, 10] and M = 20 in the region r ^ 10, 
T = 1.025 10 11 . The time-step was At = 0.01. In the top 
panel we show the results for the autocorrelation function. 
The solid (red) line shows the theoretical prediction obtained 
from Eq. (J8), while the open circles show the result of the 
numerical simulations. By performing a nonlinear fit (dashed 
blue line), the autocorrelation function shows an asymptotic 
decay compatible with a power-law r~' 3c , with (3 e — 0.79. In 
the inset of the top panel we show the numerical simulation 
(circles) relative to the mean square displacement (| A.s(£) | 2 ) . 
A nonlinear fit (solid blue line) shows that (|As(t)| 2 ) oc t s 
with 8 — 1.26. The bottom panel shows the stationary pdf of 
the process. Again the solid (red) line shows the theoretical 
prediction of Eq. (|13[1 , while the open circles show the result 
of the numerical simulations. 



tionary pdf with tails that can be gaussian and exponen- 
tial. The processes are obtained by simply performing 
a coordinate transformation of the additive process de- 
scribed in Eq. ©. Starting from such specific process, 
we have given analytical evidence that the considered 
processes have the wanted stationary pdf and we have 
given numerical evidence that the autocorrelation func- 
tion shows a power-law decay. Specifically, we find that 
for large values of time lag the autocorrelation function 
decays like (y(t)y(t + r)) oc t~@ where j3 grows with a 
with a law which seems compatible with = a/2 — r\ 
where 77 is a parameter which depends from the specific 
tails of the stationary pdf. However, when the tails are 
gaussian rj does not show any dependance from the vari- 
ance of the pdf. The above linear law holds true also in 
the case of the additive process of Eq. ([5]), with rj = 3/2. 

It is worth remarking that in principle more general 
processes can be obtained (i) by choosing different coor- 
dinate trasformations or (ii) by appropriately engineering 
the shape of the quantum potential Vs {x) of Eq. (JSJ) in 
the region [-L, L], This would result in a different shape 
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of the stationary pdf in that region. When doing that, 
the asymptotic power-law behaviour of the autocorrela- 
tion function is not modified. In this paper we preferred 
to consider a linear transformation and Vs(x) = — Vq in 
the region [— L, L] only because this allows us to analyt- 
ically obtain the eigenfunctions on the whole real axis 
and to obtain a numerical theoretical prediction for the 
autocorrelation function of the stochastic processes con- 
sidered. 

Starting from the process of Eq. ([5]), stationary pdfs 
with tails different from exponential or gaussian ones can 
be obtained by introducing appropriate coordinate trans- 
formations. In all cases the autocorrelation functions can 
be obtained, at least numerically, by using the same ap- 
proach illustrated in this paper. 

To our knowledge, this is the first evidence of power- 
law correlated stationary Markovian processes with gaus- 
sian or exponential tails in the stationary pdf. It is worth 
remarking that the existence of power-law correlated pro- 
cesses wi th g aussian tails does not contraddict the Doob 
theorem j34(, because the Doob theorem deals with the 
case when the process admits a stationary pdf and a 2- 



point conditional transition probability which are both 
gaussian on the whole real axis and non singular. In our 
case we only have gaussian tails in the stationary pdf. 

Our results help in clarifying that even in the con- 
text of Markovian processes long-range dependencies are 
not necessarily associated to the occurrence of extreme 
events. It is worth mentioning that the processes intro- 
duced in section IIIII and section IIVI are in the basin of 
attraction of the Gumbel distribution [35|, although the 
one of Eq. ([5]) is in the basin of attraction of the Frechet 
distribution. 

Moreover, our results can be relevant in the modeling 
of complex systems with long memory. In fact, processes 
with long-range interactions are often modeled by means 
of the Fractional Brownian motion (FBm), multifractal 
processes, memory kernels and other. Here we provide 
simple processes associated to Langevin equations thus 
showing that memory effects can still be modeled in the 
context of continuous time stationary Markovian pro- 
cesses, i.e. even assuming the validity of the Chapman- 
Kolmogorov equation. 
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